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Abstract 

Using  infrared  detectors,  Zehnder  and  Rosakis, [1]  Zehnder  and  Kallivayalil[2]  and  Mason 
and  Rosakis,[3]  have  recorded  the  temperature  field  around  a  dynamically  propagating  crack 
tip  in  several  metals.  At  the  same  time,  Tzou[4,5]  has  suggested  that  the  temperature  field 
around  a  propagating  crack  tip  might  exhibit  some  of  the  characteristics  of  hyperbolic  heat 
conduction.  In  this  paper  it  is  shown,  by  using  a  corrected  solution  of  the  hyperbolic  heat 
conduction  equation  for  a  traveling  point  source  and  by  using  an  experimental  estimate  of 
the  active  plastic  zone  (heat  generating  zone)  at  a  crack  tip,  that  the  effects  of  hyperbolic  heat 
conduction  around  a  propagating  crack  tip  are  not  observed  for  the  conditions  of  the  experiments. 
It  is  shown  that  due  to  the  adiabatic  conditions  at  the  crack  tip  the  solution  of  the  hyperbolic 
heat  equation  is  indistinguishable  from  the  solution  of  the  parabolic  heat  conduction  equation 
for  crack  propagation  in  steel.  Some  features  of  the  solution  are  examined  for  conditions  other 
than  those  in  the  experiments. 
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1.  Introduction 


It  is  understood  that  the  classical,  or  parabolic,  heat  conduction  equation  has  an  inherent 
pathology.  That  is,  when  a  point  source  is  introduced  into  a  conducting  medium,  the  parabolic 
heat  conduction  equation  predicts  that  its  presence  is  instantaneously  felt  throughout  the  medium. 
Often  this  pathology  is  referred  to  as  “the  infinite  speed  of  heat  propagation,”  and  it  has  been 
addressed[6]  through  the  introduction  of  a  new  heat  flux  law.  Usually,  heat  flux,  q,  is  related 
to  the  gradient  of  temperature  by  Fourier’s  Law, 

q  =  -kVT. 

When  this  relation  is  combined  with  the  expression  for  the  continuity  of  energy. 


-V  •  q  +  S  =  ()cpf . 


where  the  dot  refers  to  time  differentiation,  the  parabolic  heat  conduction  equation  results, 

aV2T-t  =  — —5  (1.1) 

pcp 

(In  this  paper  temperature,  T,  is  implicitly  taken  as  the  change  in  temperature  above  ambient, 
T  =  Tactual  —  Tmmbient  )  If,  instead  of  Fourier’s  Law,  a  new  heat  flux  law  is  used, 

^q  +  q  =  -frVT, 

then  through  similar  manipulations  the  “hyperbolic”  heat  conduction  equation  results, 

oV2T-^f~T=  — —  (S+^S)-  (1.2) 

vl  pcp  \  V*  / 

The  introduction  of  the  new  heat  conduction  law  has  been  justified  by  several  authors,  (for  a 
review  see  Tzou[4])  and  some  experimental  work  has  shown  that  this  different  heat  conduction 
equation  is  more  appropriate  at  large  distances  from  a  point  source  or  at  very  short  times  after 
the  introduction  of  a  point  source.[7]  An  estimate  for  the  speed  of  heat  propagation  in  materials, 
t>,  may  be  found, [8]  and  for  steel  this  estimate  predicts  a  speed  of  heat  propagation  on  the  order 
of  103  m/s. 


Zehnder  and  Rosakis[l]  and  Mason  and  Rosa'  is[3]  have  measured  the  temperature  at  the 
tip  of  a  dynamically  propagating  crack  for  crack  tip  velocities  ranging  from  600  to  900  m/s 
in  4340  steel.  The  temperature  fields  show  a  strong  difference  in  their  geometric  nature,  see 
Figure  1,  and,  since  these  velocities  are  comparable  to  the  estimate  for  the  wave  speed  of  heat 
propagation,  it  is  suggested  that  this  difference  may  be  due  to  hyperbolic  heat  conduction.  In 
order  to  determine  whether  hyperbolic  heat  conduction  is  active  in  these  experiments,  it  is  the 
purpose  of  this  investigation  to  produce  theoretical  temperature  fields  for  the  problem  of  a  crack 
propagating  in  a  hyperbolically  conducting  metal.  It  is  fundamentally  important  to  note  that  the 
propagation  of  a  crack  produces  plastic  deformation  at  the  crack  tip  and  that  the  plastic  work 
generated  during  this  deformation  is  mostly  transformed  into  heat.  Therefore,  the  problem  of 
calculating  the  temperature  field  around  a  propagating  crack  tip  is  more  correctly  stated  as  the 
problem  of  calculating  the  temperature  field  around  a  propagating  heat  source  zone  with  careful 
attention  paid  to  the  boundary  conditions  at  the  crack  faces. 

Several  investigations  of  the  temperature  field  around  a  propagating  crack  have  been  re¬ 
ported  for  parabolic  heat  conduction.  (For  a  complete  review  of  this  field  of  work  see  Mason  and 
Rosakis[3].)  The  most  pertinent  here  is  the  work  of  Weichert  and  Schonert.[9]  These  authors  nu¬ 
merically  calculated  the  temperature  fields  around  propagating  heat  source  zones  by  integrating 
the  solution  of  the  parabolic  heat  conduction  equation  for  a  travelling  point  source  over  a  rect¬ 
angular  heat  source  zone  of  constant  magnitude.  A  similar  methodology  is  used  here.  It  should 
be  noted  that  in  all  theoretical  investigations  of  the  temperature  field  around  a  propagating  crack 
tip  the  assumptions  made  about  the  plastic  work  zone  and,  consequently,  the  heat  source  zone 
have  proven  to  be  extremely  important  in  the  calculation  of  maximum  temperature.  Due  to  the 
complexity  of  the  three  dimensional  deformation  where  temperature  measurements  are  performed 
(at  the  surface  of  the  specimen),  it  is  difficult  to  accurately  describe  the  plastic  work  zone  in 
closed  form  and  most  assumptions  result  in  an  oversimplification  of  the  experimental  problem. 
Consequently,  in  this  work  the  shape  of  the  plastic  work  zone  is  estimated  from  the  experiments 
of  Zehnder  and  Rosakis[l]  and  Mason  and  Rosakis|3]  directly,  eliminating  over-simplification 
and,  thus,  more  closely  approximating  the  experimental  data. 


2.  Theoretical  Development 

First  it  is  necessary  to  find  the  solution  of  the  hyperbolic  heat  conduction  equation  for  a 


point  source  travelling  at  constant  velocity.  Although  earlier  attempts  to  provide  solutions  of  Eq 
(1.2)  for  a  traveling  point  source  exist,[10,l  1]  the  solution  given  here  includes  some  corrections 
that  significantly  alter  the  behavior  of  the  solution. 

We  are  interested  in  solving  the  equation, 

V\T  - —  [S(r,  -  at.:2)  +  ~  ^-2)]  ,  (2.1) 

Vz  Q  pCpO  L  V*  J 

for  a  point  source  travelling  at  a  constant  velocity,  a,  in  the  zx  direction,  i.e. 

S(z\  -  at,z2)  =  6(z i  -  df)<5(;2)- 


where 


a 


which  when  inserted  in  equation  (2.3)  yields 


v?/  -  *2/  =  -- i-exp [«£,]  { «(£,  v/TTmT)^ )  -  ■  T°W’,8-3‘  — <(6)1  ■ 

PCpQ  ^  ay/ 1  -  M2  ofi  J 

(2.5) 

The  left  hand  side  of  Eq  (2.5)  is  the  modified  Helmholtz  equation.  The  Green’s  function  solution 
of  this  equation  is  given  by  Arfken,[12] 

G(r«,r<)*-  -^A'„(K|r<-r<|)  (2.6) 

where  K0(w)  is  the  zeroth  order  modified  Bessel’s  function  of  the  second  kind  and 

|r€  -  r< |  =  \/(£i  -  Ci)'2  +  (£2  -  Cj)2- 


The  first  vector,  r*  refers  to  the  point  of  interest  while  the  second  vector,  r**,  refers  to  the 
location  of  the  point  source.  Using  the  Green’s  function  for  the  modified  Helmholtz  equation, 
the  solution  to  equation  (2.5)  may  be  found. 


/(Ci, £2)  =  7;  ■  ■  f  A’o(K|r{  -  rc|)exp(«Cij(f(C,  Vl  ~  M2 )M&)d(id(2 

Znpcpa  JA( 


1  a  A/2  f  ,  . 

-  ~ - .  rrm’»r?  /  A«(k  r*  -  r'*  )cxp[h-Ci] 

2npcpa  a VI  -  A/2  JA< 


d6(  Ci  y/\  -  M2) 

dC\ 


^CzMCi^Cz. 

(2.7) 


Equation  (2.7)  may  be  divided  into  two  integrals 

/(£. ,  £2 )  =  —  (.7.  (Ci .  £2 )  4-  Mi  1 ,  i2 ))  (2.8) 

2np(vc\ 

The  first  of  these  integrals  is  found  using  the  fundamental  property  of  the  Dirac  delta  function, 


*7i(£i  ,(2)  = 


and  the  second  may  be  found  using  the  relation 


Ao(K|r<|) 

y/l-M2' 


(2.9) 


i.e. 


*«*•«»>  -  2(i -ZWffi  {*^>  +  ^(^1)} 


(2.10) 


Substituting  Eqs  (2.9)  and  (2.10)  into  Eq  (2.8)  gives  the  result 

'«>•«*>  “  4^(1  -W»  I'"  -  W!)A’“U|r£|)  +  "2FlA''(K|ri|)}  CH) 

Remembering  our  intial  guess  for  temperature  in  Eq  (2.4),  we  arrive  at  the  solution  to  Eq  (2.3), 

T,U,.U)=  hpc  ■—  ^2j3/iexp  [-!■•{,)  {(2  -  M2)A-o(«|r'|)  +  M2^-A-1(K|r«|)}  , 

'  !  (2.12) 


where 

!r*l  =  y/fi  +  C  l 

This  is  the  solution  for  hyperbolic  two  dimensional  heat  flow  around  a  point  source  traveling 
in  a  straight  line  at  a  constant  speed,  a  <  v,  in  an  infinite  body  with  T  =  0  as  r  — ►  oo. 
As  M 2  — ♦  0  this  function  gives  the  solution  of  Carslaw  and  Jaeger[13]  for  a  traveling  point 
source  in  a  parabolically  conducting  material  under  identical  boundary  conditions.  This  solution 
differs  from  the  solution  provided  by  Tzou  (10,1 1 1  due  to  the  interpretation  of  the  argument  of 
the  modified  Bessel’s  function  in  Eq  (2.6).  In  the  work  by  Tzou,  it  appears  that  the  argument 
of  the  modified  Bessel’s  function  had  been  interpreted  as  the  difference  of  the  magnitudes  of 
two  vectors  where  it  is  actually  the  magnitude  of  the  difference  of  two  vectors.  The  additional 
multiplier  of  the  second  term,  £i/|r€|,  provides  desirable  behavior  for  the  point  source  solution 
not  found  in  the  expressions  provided  by  Tzou.  Note  that  using  the  transformation  of  coordinates 
in  Eq  (2.1)  yields 


V|T  + 


a  dT 
as/ 1  -  Ml 


1 

pc,, a 


S(£,  x/l- 3/^2  )- 


oil/2  dS(fryT-M2,fr) 
a/f-  M*  dfr 


(2. 


3) 


and  that  if  the  slope  of  the  source  zone,  S(£i  \/l  -  M1,  £2  )* is  negative  in  the  direction  there 
is  a  heating  contribution  due  to  the  second  forcing  term;  if  the  slope  of  the  source  zone  in  the 
£i  direction  is  positive,  there  is  a  cooling  contribution.  The  slope  of  a  Dirac  delta  function  in 
the  direction  is  negative  infinite  for  positive  and  positive  infinite  for  less  that  zero.[14] 
Consequently,  it  is  expected  that  there  should  be  a  heating  contribution  ahead  of  the  point  source 
that  transforms  to  a  cooling  contribution  behind  the  point  source.  The  £i/|r*j  term  changes  the 
sign  of  the  second  term  in  the  solution  to  produce  this  behavior. 

On  the  crack  faces  it  is  expected  that  the  heat  flux  out  of  the  material  is  zero,  that  is, 

_  n  f  for  =  0+  and  <  0; 

92  ~  (  for  £2  =  0”  and  <  0. 


(2.14) 


Naturally,  this  is  an  idealization  of  the  actual  boundary  condition  which  would  involve  radiation 
transfer,  convective  transfer  etc.  However,  it  is  clear  that  a  fixed  temperature  over  the  entire 
crack  faces  is  not  an  appropriate  boundary  condition  and  is  not  included  in  this  analysis.  In 
parabolic  heat  conduction,  symmetry  of  the  temperature  solution  about  the  xj  axis  is  sufficient 
to  satisfy  the  boundary  condition  at  the  crack  faces  given  by  Eq  (2.14), 


=  92(6,0)  =  0. 


However,  for  hyperbolic  heat  conduction  the  symmetry  results  in  a  first  order  differential  equa¬ 
tion,  namely. 


—ita  d(jo 

rVl  -M2dTi 


+  92(6 , 0)  =  0 


It  is  obvious  from  Eq  (2.14)  that  q2  =  0  is  the  solution  to  this  equation.  However,  it  is  also 
clear  that  the  possibility  exists  for  other  flow  conditions  to  be  prescribed  on  the  crack  faces  if 
boundary  conditions  other  than  Eq  (2.14)  exist.  In  hyperbolic  heat  conduction  symmetry  about 
the  6  axes  does  not,  by  itself,  imply  that  q,  is  identically  zero  on  the  axis. 


2.2  For  the  case  M2  >  1 

By  employing  the  transformation. 


c  _  -1  -  at 

6  —  /TTT 


y/AP- T, 

6  =  -2 


Eq  (2.2)  may  be  reformulated  as 

d2T  d2T  a  dT 
+  —  + 


1 


pcpo 


a  v  A/2  -  1  ail 


We  guess  a  solution  to  Eq  (2.15)  of  the  form 

Tpttuh)  =  exp  [^6  )/(6>  6) 


(2.15) 


(2.16) 


where 


and  substitute  above  resulting  in 


f/_^Z_Kv  = 
atf  an  } 


1  f  /— -  nil/2  dSihVM2  -  1) 

-exp [-«{,]  MC,  y/M  -DM6)-  ■  A/m  36 


(2.17) 


M6) 


pCpQ 


The  left  hand  side  of  this  equation  is  the  telegraph  equation.  The  Green’s  function  for  this 
equation  may  be  found  in  Duff  and  Naylor)  15] 


G(r«,rc)  =  ij0(K|rc  -  rc|)ff(|r<  -  rc| 2)H(Ci  -  6) 


(2.18) 


where  Jo  is  the  zeroth  order  modified  Bessel  function  of  the  first  kind,  H(x)  is  the  Heaviside 
function  and 

|r*  -  r*’|  =  \/(£i  ~  Ci ~  ( £2  -  C2  )2 • 


Because  of  the  Heaviside  terms,  the  solution  in  Eq  (2.18)  is  zero  outside  a  triangular  regime 
behind  the  point  source.  There  exists  another  Green’s  function  that  is  non-zero  ahead  of  the 
point  source,  but  this  function  is  not  used  due  to  obvious  physical  arguments.  Because  the  point 
source  is  travelling  faster  than  the  wave  speed  of  heat  propagation  in  the  material,  it  is  physically 
impossible  to  have  a  finite  temperature  ahead  of  the  point  source.  Also,  due  to  the  Heaviside 
terms,  the  temperature  exhibits  a  jump  or  shock  along  lines  inclined  at  an  angle  u  to  the  negative 
£1  axis.  The  shock  angle,  u\  is  given  by  the  familiar  formula. 


u;  =  ±sin  1(1/.AJ). 


(2.19) 


Not  surprisingly,  these  thermal  shocks  are  quite  similar  to  the  shocks  found  in  super-sonic  fluid 
flow.  This  is  to  be  expected  since  Eq  (1.2)  is,  in  fact,  the  damped  wave  equation. 

Convolution  of  the  Green’s  function  for  the  telegraph  equation  with  the  forcing  term  in 
Eq  (2.15)  results  in  two  integrals.  As  before,  the  first  of  these  integrals  is  found  using  the 
fundamental  property  of  the  Dirac  delta  function. 

The  second  integral  is  found  using  the  additional  relation  6(x  -  r)S(x)  =  0, 

*«..&>  ■  -wJTZwi  ^/,<^l)}g(|r<ng(-{,). 


Combining  these  two  gives  the  result 


/«>.&>  =  4pCpa(ip-1)3/1  {(W2  -  tf(|r'|2  )*(-£,)■ 

Remembering  our  intial  guess  for  temperature  in  Eq  (2.16),  we  arrive  at  the  final  result, 

7VUi,6)  =- - — 77572  exP  Kil 

4pcpa(M2  —  l)3/2 

|(M2  -  2)/0(K|r«|)  -  Mi^/1(K|r«|)|  tf(|r<|2  )#(-£, ). 


This  is  the  solution  of  Eq  (2.3)  for  a  point  source  travelling  at  constant  velocity  with  M2  >  1. 
It  differs  significantly  from  the  solution  reported  by  TzouflOJ  1]  due  to  the  use  of  the  correct 
Green’s  function  solution  for  the  telegraph  equation.  It  must  be  re-emphasized  that  the  definition 
of  the  argument  of  the  modified  Bessel’s  functions,  |r*|,  is  significantly  different  in  this  case  than 
in  the  previous  case,  and,  furthermore,  that  the  sign  of  the  argument  of  the  leading  exponential 
term  has  changed.  Although  the  two  solutions  appear  similar,  they  are,  in  fact,  not  very  similar 
at  all;  for  example,  this  solution  does  not  converge  to  the  solution  of  Carslaw  and  Jaeger[13] 
under  any  circumstances.  Note  the  effect  of  the  cooling  term  in  the  forcing  function  on  the 
right  hand  side  of  equation  (2.15),  when  M1  <  2  the  temperature  at  the  location  of  the  point 
source  is  negative  due  to  the  fact  that  J|(tr)  — »  0  as  ir  — ►  0  while  Iq{w)  -*  1  as  w  — »  0.  If 
the  cooling  term  in  the  forcing  function  of  Eq  (2.15)  were  neglected  this  effect  would  not  be 
observed.  When  M2  >  2,  the  solution  is  always  positive. 


23  Integration  over  the  Source  Zone 

For  the  sake  of  simplicity  the  heat  source  zone,  S(£ ,  .£2 ).  has  been  assumed  to  be  defined 
only  on  the  rectangle  £1  €  (0, 6/V 1  -  M2}  and  £2  €  [-ft,  f6]  where  /  is  the  aspect  ratio  of  the 
zone;  elsewhere  it  is  identically  zero.  From  the  experimental  work  of  Zehnder  and  Rosakis[l] 
and  Mason  and  Rosakis[3]  it  is  seen  that  a  close  approximation  to  the  experimental  heat  source 
zone  is  given  by 

1  +  cos0f).' 


Sii  1,6)  J_ 

Q  2/ 


1  —  cos 


2jt>/1  - 


where  /  =  0.05.  A  comparison  of  the  approximate  plastic  zone,  S(£  1,^2).  with  the  experimen¬ 
tally  estimated  plastic  zone  is  shown  in  Figure  2  for  both  testing  conditions  .  The  inclusion 
of  1/2/  in  the  denominator  gives  this  function  characteristics  of  a  Dirac  delta  function  in  the 
£2  (or  X2»  given  below)  direction  as  /  — +  0.  Consequently,  this  function  is  also  normalized  so 
that  its  integral  in  two  dimensions  is  unity  giving  usable  results  if  /  is  allowed  to  be  zero.  The 
multiplier,  Q,  is  found  experimentally  from  the  maximum  value  of  the  work  rate  density  and 
the  relation  between  Q  and  the  maximum  of  the  function. 


Q  = 


Once  this  relation  for  the  heat  source  zone  has  been  assumed,  letting 

77  £1 


•*'  1  =  Ci 

0 

.  £.1 

.To  =  (,  >  —  —  . 

*  d 


(2.21) 


the  temperature  field  due  to  the  heat  source  zone  ahead  of  a  crack  tip  may  be  found  by  a 
convolution; 


T(; rj.x-2)  =  Qt> 


k2 


rl  rf  SlCuOi) 


a 


Q 


Tn(  X 1  -  Cl  *  -J‘2  -  C2  )rfC2^Cl  1 


or,  more  specifically,  for  M 2  <  1 

Q62 


T(x !,X2) 


47rpcro(l  —  A/2)3/2 


/'  5(c..c2)...„ 

'^s 

1 

1 

« _ 

1  ' 

O 

> 

1  -  A/2  J 

(2.22a) 


where 


a  6 


V'  =■ 


2o 


(2.22b) 


lrx  -  rM 


=  n/(.t,  -Ci)2  +  (1-JI/2)(x2-C2)2 

Normalizing  temperature  with  respect  to  Qbjpcpit  for  numerical  integration  yields 


0(xi,*2)  PqsT(Xi'X2)  L  Lf 


S(  Ci.Cz) 

Q 


exp 


-rjjn  -  Ci  )1 

1  -  A/2 


(2.23a) 


Similarly  for  M2  >  1 


6(xux2)  = 


2(A/2 

|(M2 


1  S(  Ci,C>) 


Q 


-exp 


-nix  1  -  G  ) 


M2  - 1 


2)/0 


c'’|rx  -  r*'! 


.l/2  -  1  )  |rx  - 

tf(|rx-r‘-‘|2)tf(  — r,  +C.)AV/Ci- 


V’  |  rx  -  r“*| 
M2  -  1 


(2.24b) 


Numerical  integration  of  the  solution  for  M2  <  1  proceeds  without  avail.  The  singularity 
in  the  modified  Bessel’s  function  of  the  zeroth  order  at  the  origin  is  logarithmic,  therefore 
the  integral  is  finite.  The  singularity  in  the  modified  Bessel’s  function  of  the  first  order  at 
the  same  point  is  of  order  l/|rx  -  r4*  |  which  is  integrable  in  two  dimensions.  (This  is  easily 
shown  by  a  conversion  to  polar  coordinates.)  Multiplication  of  A'i(V’|rx  —  r*»|)  by  the  factor 
(ii  —  (i)/|rx  -  r^l  does  not  change  the  order  of  the  singularity  of  this  term  since  at  ( x\  ~Ci)  = 
(x2  —  C2)  =  0  this  term  is  finite.  For  M2  >  1,  the  integration  is  quite  simple.  Since  both  modified 
Bessel  functions  of  the  first  kind  are  finite  at  the  origin,  the  only  singularity  occurs  in  the  second 
term  and  this  remains  integrable.  (This  is  easily  shown  using  polar  coordinates.) 

The  convolution,  Eq  (2.22a),  is  evaluated  numerically  for  li-  6  {.01,1,100]  and  M2  6 
[0,  .5,  .9, 100, 104].  The  modified  Bessel’s  functions  are  evaluated  either  by  using  IMSL  sub¬ 
routines  for  small  arguments  or  by  using  the  asymptotic  expressions  given  below  for  larger 
arguments.  To  find  the  value  of  the  integral  a  two  dimensional  Gauss-Legendre  scheme  is 
employed.  The  number  of  integration  points  increased  with  the  value  of  r/’* 


2.4 


Asymptotic  Analysis  of  the  Integration 
By  employing  the  transform 

r,  -  cit 


in  Eqs  (1.1)  and  (1.2)  the  following  expressions  are  found; 


J-v^r+2 1 

2>'-  *  dr, 


opcp 


and 


1  ,  M2  d2T  dT_  = 

2 1  *  2c  dr2  + 


(2.25a) 


Of)Cp 


(2.25b) 


respectively.  Note  that  as  the  parameter  y'  gets  large  the  left  hand  sides  of  both  equations  become 
more  adiabatic  and  that  the  solution  to  the  adiabatic  equation  is  given  simply  as; 


0a(*l,*2) 


ij  { 1  +  co^TT^f  ))  for  .rj  <  0  and  x2  6  [-/,/], 

'  2/  (1  +cos(7rif ))  (!  _  1  i  +  y;sin(27r.ri ))  for  aj  6  [0,1],  x2  €  [-/,/], 


v  0 


otherwise. 


(2.26) 


where  the  temperature  has  been  normalized  by  Q6/pc}>c.  From  the  definition  of  xp,  it  is  seen  that 
one  of  three  things  may  be  occurring  when  ir  — ♦  oc:  n  — *  0,  a  — *  oo  or  6  — ►  oo.  In  Eqs  (2.25a) 
and  (2.25b)  as  yb  —*  oo  the  left  hand  side  becomes  the  adiabatic  equation,  but,  depending  upon 
how  rp  is  increasing,  through  o,  a  or  <^,  the  right  hand  side  may  be  affected  as  well.  Normalizing 
the  the  temperature  with  respect  to  Q6fpcva  in  Eq.  (2.24b)  removes  this  ambiguity  in  the  limit, 
and,  thus,  the  point  source  term  in  the  integrand  of  Eq  (2.24b)  rather  than  Eq  (2.22a)  is  examined 
asymptotically  for  M2  <  1. 

For  large  arguments  the  modified  Bessel’s  functions  can  be  approximated  by 
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(2.27) 


Using  only  the  first  terms,  the  point  source  in  the  integrand  for  M2  <  1,  Eq  (2.24b),  can  be 
rewritten  as 
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where 

|rx|  =  y.rf  +  d-^J.r2. 


Given  that  for  small  .r2/.r1 
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(2.29) 


Eq  (2.28)  may  be  expanded  to  take  the  following  form. 


2sgn(.ri  )A/2(  1  -  M2 )  +  (2  -  M-  +  sgn(.r i  )M2)(  1  —  A/2)  / £2  \2 

4  V 


When  x\  >  0  the  exponential  term  exp  [— i,’(  ri  +  |.i'i  | )/( 1  -  M2)]  dominates  the  solution  and 
Tp(x \,x2)  ~  0;  on  the  other  hand,  when  ,r,  <  0,  (.i-j  +  |.-r  1 1)  =  0,  the  same  exponential  term 
disappears  and  sgn(.rj )  is  always  negative.  Thus, 


0P(xi,x2) 


for  a-j  <  0; 
for  x\  >  0 


(2.30) 


It  is  clear  that,  as  i/’/2\xl\  becomes  large,  the  second  term  in  the  expansion  becomes  negligible 
and  the  solution  becomes  independent  of  M2  Consequently,  the  solution  for  the  hyperbolic  heat 
conduction  equation  and  the  solution  for  the  parabolic  heat  conduction  equation  converge  to  the 
same  result.  It  is  also  clear  that  as  ?/’/2|.r  1 1  — ♦  oc,  0,,(.r].: r2)  becomes,  as  expected  from  the 
adiabatic  solution,  a  Dirac  delta  function|12|  in  v-2.  However,  contrary  to  the  solution  for  the 
adiabatic  equation  and  in  keeping  with  the  boundary  condition  T  — ►  00  as  rx  — »  00,  the  function 
in  Eq  (2.30)  loses  its  Dirac  delta  function  character  as  |.T]  |  grows. 

Not  surprisingly,  similar  results  are  found  for  the  case  M2  >  1.  Asymptotic  analysis  using 
the  approximations 


and  following  the  same  steps  outlined  above  leads  to  the  same  leading  term  characterization  of 
the  solution  as  for  the  case  M2  <  1.  Since  the  governing  equation  becomes  adiabatic  as  the 
approximation  becomes  more  accurate  it  is  to  be  expected  that  both  solutions  converge  to  the 
same  result  as  they  do. 


3.  Results  and  Discussion 

The  solution  for  the  simplest  case,  the  adiabatic  case  in  Eq  (2.26),  is  shown  in  Figure 
3.  Behind  the  heat  source  zone  the  adiabatic  solution  is  characterized  by  contours  of  constant 
temperature  extending  from  the  heat  source  zone  to  .i  j  =  —  oo  parallel  to  the  crack  faces.  In 
front  of  the  crack  tip,  in  the  heat  source  zone,  the  gradient  of  temperature  along  lines  extending 
radially  from  the  crack  tip  appears  to  be  constant.  Although  it  isn’t  true  mathematically  that  the 
radial  gradient  in  the  heat  source  zone  is  constant,  it  is  noted  that  the  expression  for  the  zone  used 
here  is  an  approximation  to  the  experimental  data  and  that  “averaging”in  the  experiments  due  to 
the  finite  infrared  detector  size  reduces  the  accuracy  of  the  measurements  near  the  crack  tip.[2] 
Generally,  it  is  expected  that  if  adiabatic  conditions  prevail  at  the  crack  tip,  a  nearly  constant 
radial  gradient  may  be  recorded  experimentally.  Good  quantitative  agreement  is  found  between 
the  predicted  temperature  rise  and  the  measured  temperature  rise.  A  maximum  temperature  rise 
of  298°C  is  seen  in  the  experiment  when  a  =6(X)  m/s,  and  the  predicted  result,  304°C  occurring 
along  the  crack  line,  is  very  close  to  that  measurement.  The  minimum  temperature  is  0°C 
occurring  everywhere  outside — except  directly  behind — the  heat  source  zone. 

The  results  of  the  numerical  integration  for  the  M2  <  1  cases  are  plotted  in  Figures  4 
through  6.  For  small  i p  it  can  be  seen  in  Figure  4  that  there  is  a  marked  dependence  of  the 
temperature  field  upon  M2.  Most  notably,  for  M  ’  e  [.5.  .9]  a  region  of  temperature  drop  is 
seen  behind  the  source  zone.  This  is  an  interesting  but  unrealistic  effect  since  the  combination 
of  parameters  in  these  plots,  r  =  .01  and  M2  £  (.5.  .9],  is  unlikely  to  be  seen  experimentally. 
As  M 2  — +  1  the  temperature  field  becomes  more  localized  near  the  origin  and  the  temperature 
at  any  given  point  behind  the  heat  source  zone  decreases  with  increasing  M2.  For  M2  =  0  the 
maximum  normalized  temperature,  6maz,  occurs  roughly  at  the  maximum  of  the  heat  source  zone 
and  is  equal  to  .1%  of  the  adiabatic  maximum.  The  minimum,  naturally,  is  zero  as  r  — ►  oo.  For 
M 2  =  .5  the  maximum  is  moved  forward  to  roughly  the  location  of  the  minimum  slope  of  the 
source  zone  where  8moz  is  equal  to  1%  of  the  adiabatic  maximum.  The  minimum  occurs  roughly 
at  the  location  of  the  maximum  slope  of  the  source  zone  where  the  normalized  temperature  is 
equal  to  -1%  of  the  maximum  adiabatic  temperature.  For  the  case  M2  =  .9  the  location  of  the 
maximum  and  minimum  are  not  changed  with  respect  to  the  case  M2  =  .5  but  the  magnitude 
does,  8maz  is  equal  to  10%  of  the  adiabatic  maximum  while  0„,,„  is  equal  to  -10%  of  the 
adiabatic  maximum.  The  behavior  seen  for  the  two  cases  when  M2  ^  0  is  exactly  as  expected 


from  Eqs  (2.13)  and  (2.12).  When  M2  ±  0  there  is  a  contribution  due  to  the  slope  of  the  heat 
source  zone  as  seen  in  Eq  (2.13).  This  contribution  has  a  cooling  effect  when  the  slope  is  a 
maximum  and  a  heat  effect  when  the  slope  is  a  minimum;  as  M 2  — *  1  the  significance  of  the 
heating  and  cooling  due  to  the  modified  Bessel’s  function  of  the  first  order  becomes  greater 
and  more  dominant.  Thus,  as  M2  — »  1,  the  maximum  temperature  moves  to  the  location  of 
minimum  slope  and  the  minimum  temperature  moves  to  the  location  of  the  maximum  slope. 
Both  maximum  and  minimum  increase  in  magnitude  as  M2  — >  1. 

For  larger  Figure  5,  the  dependence  of  the  temperature  field  upon  M2  begins  to  disappear. 
No  negative  temperature  changes  are  seen,  however,  when  M2  =  .9  the  temperature  does  exhibit 
a  positive  minimum  near  the  tail  end  of  the  heat  source  zone.  Thus,  some  cooling  effects  due  to 
the  second  source  term  in  Eq  (2.13)  remain.  Although,  far  from  the  heat  source  zone,  x\/b  >  1, 
all  three  fields  are  equal.  As  M2  ->  1  a  localization  of  the  temperature  near  the  source  zone 
is  still  seen.  The  maximum  temperature  is  located  near  the  maximum  of  the  source  function  in 
each  sub-case  although  it  moves  forward  slightly  as  the  thermal  mach  number  increases.  The 
maximum  normalized  temperature  increases  with  thermal  mach  number;  8maz  =  .75,1,1.5  for 
M2  —  0,  .5,  .9  respectively.  These  values  are  3.25%,  5%  and  7.5%  of  the  maximum  temperature 
under  adiabatic  conditions. 

For  even  larger  xj\  Figure  6,  the  dependence  upon  M2  disappears  completely  as  expected 
from  the  asymptotic  analysis,  Eq  (2.30),  and  the  hyperbolic  heat  conduction  solution  is  indis¬ 
tinguishable  from  the  solution  for  parabolic  heat  conduction  at  all  values  of  M2  <  1.  No 
localization  of  the  temperature  near  the  heat  source  zone  is  discemable.  The  maximum  temper¬ 
ature  occurs  slightly  behind  the  maximum  of  the  heat  source  zone  with  value  of  7.5 — 32.5%  of 
the  maximum  temperature  under  adiabatic  conditions. 

Small  V’  accounts  for,  a  small  source  zone,  6,  low  velocity,  a,  or  a  large  thermal  diffusivity, 
a.  (See  Eq  (2.22b).)  For  metals  these  conditions  do  not  reflect  the  usual  experimental  conditions. 
For  4340  steel,  a  =  10-5,  and  af>  2  for  both  experiments  shown  in  Figure  1  giving  r />  ~  105  in 
both  cases.  For  comparison  the  results  of  the  integration  for  =  105  have  been  plotted  in  Figure 
7.  It  is  seen  that  the  theoretical  temperature  field  matches  the  experimental  results.  Figure  1,  well 
for  a  =  600  m/s  and  that  both  the  theoretical  temperature  field  and  the  experimental  temperature 
field  resemble  the  adiabatic  solution.  Figure  3.  Contours  extend  from  the  heat  source  zone  nearly 
parallel  to  the  crack  faces  toward  .tj  =  -oc  before  curving  in  to  meet  the  negative  xi  axis. 


The  results  for  xl>  =  105  differ  from  the  adiabatic  solution  in  that  the  temperature  decreases  as 
|si  |  increases  behind  the  heat  source  zone.  It  is  reiterated  that  this  decrease  is  expected  due  to 
the  zero  temperature  boundary  condition  at  |.ri|  =  oc.  This  boundary  condition  more  closely 
approximates  the  experimental  condition  than  an  adiabatic  boundary  condition  because  there  is, 
after  all,  a  small  but  significant  amount  of  heat  conduction  leading  to  a  cooling  of  the  specimen 
as  |xi  |  — »  oo.  It  is  noted  that  the  temperature  decreases  more  rapidly  behind  the  crack  tip  in 
the  experiments  than  in  the  theory.  This  is  attributed  to  heat  loss  by  radiation  and  convection 
at  the  surface.  Kuang  and  Alturif  3  6]  have  included  these  effects  in  there  numerical  parabolic 
heat  conduction  analysis,  and  they  report  a  more  rapid  decrease  in  the  temperature  behind  the 
crack  than  when  radiation  and  convection  are  not  neglected.  For  a  =  900  m/s  in  Figure  1  a 
discrepancy  with  the  theoretical  results  for  M-  <  1  in  Figure  7  is  seen.  The  contours  seen 
experimentally  emanate  from  the  source  zone  outwardly  before  curving  in  to  meet  the  negative 
Xj  axis.  (The  predicted  maximum  temperature  at  the  crack  tip  is  in  good  agreement  with  the 
experimental  measurement,  however,  owing  to  the  adiabatic  conditions  at  the  crack  tip.) 

For  M2  >  1  an  exemplary  plot  is  shown  in  Figure  8.  The  resemblance  of  this  figure  to 
the  results  recorded  in  Fig  1  for  a  =  900  m/s  is  striking,  however,  the  experimental  temperature 
field  for  a  =  900  m/s  shows  what  might  be  interpreted  to  be  thermal  shocks  at  a  shock  angle 
corresponding  to  a  large  mach  number,  M  100.  If  this  value  for  M  is  accurate  then  when 
a  —  600  m/s  the  temperature  field  should  also  show  similar  thermal  shocks;  it  doesn’t.  Also, 
because  r/’  =  105  experimentally  not  103  as  in  Fig  8,  the  theoretical  temperature  field  for 
thermally  super-sonic  crack  propagation  at  a  =  900  m/s  actually  resembles  the  temperature  field 
in  Fig  7.  This  is  expected  from  the  asymptotic  analysis  and  has  been  checked  numerically  by 
the  authors.  Therefore,  it  is  clear  that  no  evidence  of  thermally  super-sonic  crack  propagation  is 
observed  in  these  experiments. 

The  difference  in  between  the  experimentally  observed  temperature  field  at  a  =900  m/s 
and  the  predicted  field  at  the  same  crack  speed  is  due  to  crack  face  opening[3].  The  opening 
velocity  required  to  produce  the  observed  effect  is  7.5  m/s,  a  reasonable  crack  opening  velocity. 
It  is  expected  that  the  crack  face  opening  speed  depends  upon  the  initial  static  stress  intensity 
factor  and  the  crack  velocity.  From  simple  theory  one  might  expect  that  the  crack  face  opening 
velocity  increases  by  as  much  as  5  times  when  the  crack  speed  increases  from  600  m/s  to  900 
m/s.[17]  This  suffices  to  explain  the  difference  in  the  temperature  fields  in  Figures  1  and  7. 


4.  Conclusions 

From  the  asymptotic  analysis  of  the  traveling  point  source  solution  of  the  hyperbolic  heat 
conduction  equation,  it  is  clear  that,  for  crack  propagation  in  metals  when  the  crack  speed  is 
either  lower  or  higher  than  the  material  heat  propagation  speed,  the  difference  between  hyperbolic 
heat  conduction  and  parabolic  heat  conduction  is  negligible  (See  Figures  6  and  Figure  7).  As  the 
factor  rp  =  a<5/2a  gets  large,  (</>  ~  105  for  a  crack  propagating  in  4340)  solutions  for  a  traveling 
point  source  in  a  hyperbolic  or  parabolic  material  converge  to  the  same  result.  This  end  result 
is  insensitive  to  changes  in  the  thermal  mach  number,  M,  and  is  very  similar  to  the  solution  for 
a  traveling  source  zone  in  an  adiabatic  material. 

The  temperature  field  in  Figure  1  exhibits  quasi-adiabatic  heat  conduction  behavior  when 
a  =  600  m/s.  Contours  of  constant  temperature  extend  from  the  crack  tip  to  |.tj|  =  — oo  nearly 
parallel  to  the  crack  faces  before  curving  in  to  meet  the  negative  rj  axis.  A  region  of  nearly 
constant  radial  temperature  gradient  is  observed  ahead  of  the  crack  tip.  The  temperature  field 
greatly  resembles  the  solution  for  adiabatic  conditions  due  to  the  high  crack  velocity,  a,  small 
heat  production  zone,  6,  and  low  thermal  diffusivity,  n,  of  4340  steel. 

When  a  =  900  m/s  the  temperature  field  does  not  show  behavior  that  is  predicted  by 
parabolic  or  hyperbolic  heat  conduction.  It  is  understood  that  this  difference  is  due  to  opening 
of  the  crack  faces  behind  the  crack  tip.[3|  Thus,  no  evidence  of  hyperbolic  heat  conduction  is 
observed,  and  it  is  concluded  that  the  traveling  point  source  shows  little  promise  as  an  experi¬ 
mental  method  for  the  investigation  of  the  hyperbolic  heat  conduction  effect.  A  feasible  source 
zone  size  combined  with  the  normal  range  of  material  parameters  for  engineering  materials  con¬ 
sistently  leads  to  near  adiabatic  conditions  around  the  source  zone  unless  a  very  low  thermal 
wave  speed  is  expected  (~  1  m/s). 
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Temperature  fields  around  a  crack  propagating  in  oil-quenched  4340 
steel  at  two  different  velocities.  The  maximum  temperature  at  higher 
velocities  is  450°C|  1  ]  while  the  maximum  temperature  at  600  m/s  is 
300°C[3]. 

Experimental  measurements  and  theoretical  approximations  to  the  plas¬ 
tic  work  zone  on  the  surface  of  oil-quenched  4340  steel  during  dy¬ 
namic  fracture.  At  900  m/s  the  work  zone  is  smaller!  1]  than  the  work 
zone  at  600  m/s.|3|  For  each  experiment  the  aspect  ratio  is  the  same, 
/  =  .05,  but  the  size,  b  differs.  At  900  m/s  b  =  5.0;  at  600  m/s 
b  =  7. 5. (See  section  2.3.) 

The  temperature  field  due  to  the  approximate  plastic  work  zone  shown 
in  Figure  2  for  adiabatic  conditions.  The  box  indicates  the  region 
of  non-zero  plastic  work.  (See  section  2.3.)  Good  qualitative  and 
quantitative  agreement  with  the  experimental  results  for  o  =  OOOm/s 
are  seen.  The  predicted  temperature  rise  for  <)  =  900 m/s  is  also  good, 
however,  the  shape  of  the  field  does  not  agree  with  that  measured 
when  <i  =  900m/.'. (See  Figure  1.) 

The  normalized  temperature  field  for  a  propagating  source  zone  with 
tl'  =  .01  and  M2  e  [0.  5.  9j.  (See  Eqs  (2.22a)  and  (2.22b).)  Note 
that  the  existence  of  a  temperature  drop  is  exhibited  for  M2  €  [.5,  .9] 
and  that  the  temperature  field  becomes  more  localized  around  the 
source  zone  with  increasing  .\/-\ 

The  normalized  temperature  field  for  a  propagating  source  zone  with 
xi'  =  1  and  M2  €  [0.  .5.  .9).  There  exist  no  temperature  drops  as  for 
r  =  .01  in  Figure  4..  however,  the  temperature  field  becomes  increas¬ 
ingly  localized  around  the  heat  source  as  the  thermal  mach  number, 
M,  increases  leading  to  higher  temperature  rises  at  the  maximum. 

The  normalized  temperature  field  for  a  propagating  source  zone  with 
xy  =  100  and  M2  €  [0.  9J.  The  two  fields  are  virtually  indistinguish¬ 
able.  No  dependence  upon  .1/-  is  seen. 

The  theoretical  temperature  fields  for  approximate  experimental  con¬ 
ditions  when  <7=  900  m/s  and  600  m/s.  Good  agreement  is  seen  in 
the  predicted  maximum  temperature  at  the  crack  tip,  however  some 
discrepancies  occur  between  the  general  shape  of  the  field  in  this  fig¬ 
ure  and  the  temperature  fields  shown  in  Figure  1.  (x!'  =  105,  6  =  7.5 
mm  for  n  =  000  m/s  and  b  =  5  mm  for  n  =  900  m/s.) 

The  theoretical  temperature  field  when  M2  =  104  and  xp  =  103  re¬ 
sembles  that  in  Fig  1  when  a  =  900  m/s.  This  resemblance  is  coin¬ 
cidental,  however,  because  the  large  value  of  M2  indicates  thermally 
supersonic  behavior  should  also  be  expected  when  a  =  600  m/s  and 
no  such  behavior  is  observed. 
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heat  flux  vector 
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